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ABSTRACT 

Topological defects are thought to be left behind by the cosmolog- 
ical phase transitions which occur as the universe expands and cools. 
Similar processes can be studied in the phase transitions which take 
place in the laboratory: "Cosmological" experiments in superfluid He^ 
and in liquid crystals were carried out within the past few years, and 
their results shed a new light on the dynamics of the defect-formation 
process. The aim of this paper is to review the key ideas behind 
this cosmology - condensed matter connection and to propose new 
experiments which could probe heretofore unaddressed aspects of the 
topological defects formation process. 



1. INTRODUCTION 

"Big Bang" leaves a radiation-dominated universe at a temperature close to the 
Planck temperature. As the resulting primordial fireball expands, temperature falls, 
thus precipitating a sequence of symmetry-breaking phase transitions, the last of which 
converts quark-gulon plasma into nucleons. In the intervening period of expansion, 
many such transitions could have taken place. The exact number as well as the nature 
of these transitions is not known, but it is often suggested that the origin of the 
observed structures — galaxies, clusters, large-scale voids — is intimately related to the phase 
transitions which may have taken place very early in the history of the universe. For 
example, inflationary cosmology (Linde, 1990) traces the origins of the universe "as we 
know it" to the era in which the very first symmetry breaking in the GUT scheme was 
taking place. Indeed, in many models inflation happens in the course of the (usually 
second-order) phase transition (Guth, 1981; Albrecht and Steinhardt, 1982; Linde, 1982), 
leaving behind density perturbations which can act as seeds of the structures seen today 
in the universe. 

In addition to inflation there may have been many other phase transitions, which 
could have left their imprints on the mass distribution of the universe. In particular, for 
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a large class of field theories the state left after the phase transition contains topological 
defects — fragments of the original high-energy pre-transition (symmetric) "false vacuum" 
which became locked in within the topologically stable configurations of the new, broken 
symmetry "true vacuum" (Zeldovich, Kobzarev, and Okun, 1974; Kibble, 1976). Such 
frozen-out relics of the distant past come in three basic varieties: monopoles, (cosmic) 
strings, and domain walls. The kind of defects which will be left behind by the transition 
depends on the interplay between the symmetry which is being broken during the phase 
transition (more specifically, on the homotopy group of the manifold of the degenerate 
broken symmetry vacua) and the symmetries of the space in which the evolution is taking 
place (Kibble, 1976). 

Monopoles are an inevitable outcome of the GUT-era phase transition and a major 
threat to cosmological models. Their mass would dominate the universe and force it to 
recoUapse well before the present era. Infiationary cosmology was in part motivated by the 
desire to dilute monopoles (Guth, 1981). Cosmic strings are "optional" and may result in 
density perturbations not too different from those produced by infiation (Vilenkin, 1985). 
Domain walls are again a threat — they are simply too massive to be consistent with either 
the observed structures or the (relatively small) perturbations of the cosmic microwave 
background. 

The process of formation of topological defects — the main subject of this paper — is 
perhaps the least understood (and most interesting) aspect of their evolution. The nature 
of this process depends on the nature of the phase transition. In the course of the first-order 
phase transitions, a new broken symmetry phase usually nucleates after the temperature 
falls some distance below the critical temperature T^. Separated islands of the new phase 
form independently and expand, resulting in a local selection of the broken symmetry 
vacuum. As a result of these independent selections, the final configuration may not be able 
to get rid of the locked-out fragments of the original, pre-transition vacuum. First-order 
phase transitions can also occur through the process known as spinodal decomposition 
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which by-passes the post-nucleation epoch of coexistence of the two phases. 

Second-order phase transitions are still more interesting. There the phase transfor- 
mation occurs (approximately) simultaneously throughout the volume. However, unless 
the critical temperature is traversed infinitesimally slowly, the resulting broken symmetry 
phase will contain many distinct regions with different choices of the local vacuum. This 
is because the selection of the new vacuum has to be somehow communicated if the same 
choice is to be made elsewhere, and the speed at which this information can be propagated 
is finite. Hence, appearance of the topological defects is — as was observed by Kibble in his 
seminal paper (Kibble, 1976) — a direct consequence of causality. 

In the cosmological setting the ultimate causal limit on the size of the approximately 
uniform patches of the new vacuum is set by the size of the horizon at the time when 
the transition is taking place (Kibble, 1976). In the laboratory this ultimate limit is not 
relevant. Indeed, also in the cosmological context the process which will set the size of 
the uniform patches of the broken symmetry vacuum is likely to be associated with the 
dynamics of the order parameter rather than with the causal horizon (which can supply but 
an ultimate upper limit). In the original paper. Kibble appealed to the process of thermal 
activation, which is thought to occur in the course of the second-order phase transitions 
below the critical temperature, Tc, but above the so-called Ginzburg temperature Tq- 
When the temperature T of the system falls in between Tq and Tg, free energy required 
to "flip" locally uniform (correlation length-sized) volume of the new phase is less than 
the available thermal energy per degree of freedom, ~ ksT. Hence, while locally the 
new vacuum has been selected, this choice is not really permanent — transitions between 
different possible vacua occur on a relaxation time scale, which can be much shorter than 
the time over which temperature traverses the region between and Tq- Consequently, 
it was originally postulated that the initial density of topological effects is set by the 
correlation length at the Ginzburg temperature, C(^g) (Kibble, 1976). 
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In the proposal for the early universe-like experiments suggested as a superfluid test 
of the cosmological scenario (Zurek, 1984; 1985; 1993) a different alternative was put 
forward. It focuses on the nonequilibrium aspect of the phase transition and predicts that 
the characteristic domain size should be set by the critical slowing down. Order parameter 
can adjust only on a relaxation time scale which diverges at T = Tc. Thus, as the critical 
temperature is approached from above, at a certain instant t evolution of the perturbations 
of the order parameter will become so sluggish that the time spent by the system in the 
vicinity of Tc will be comparable with the relaxation time scale itself. The correlation 
length corresponding to such a "freezeout instant" will set the size of the regions over 
which the same vacuum can be selected. Hence, it will set the resulting density of the 
topological defects. This prediction (Zurek, 1984; 1985; 1993) has been now borne out by 
the experiment (Hendry et at, 1994; Hendry et al., 1995; Zurek, 1994). 

The aim of this paper is to discuss the dynamics of the non-conserved order parameter 
during the course of the phase transitions in the wake of the "cosmological" quench 
experiments which were recently carried out in superfluid helium (Hendry et ai, 1994; 
Hendry et al., 1995) and in liquid crystals (Chuang et al., 1991; Yurke, 1995; Bowick et al, 
1994) and to propose further experiments which can shed a new light on the heretofore 
unaddressed aspects of the symmetry breaking. We shall focus on second-order phase 
transitions in superfluids and in superconductors, which constitute two experimentally 
attractive alternative implementations of symmetry breaking phase transitions in systems 
with a global and local gauge, respectively. 

In the next section we shall briefly review "static" (equilibrium) aspects of symmetry 
breaking in field theories, in the superfiuid helium, and in the superconductors. Section 3 
will discuss topological defect formation in superfluid He^, compare theoretical predictions 
with the recent experiment, and consider prospects for analogous experiments in supercon- 
ductors. Section 4 will focus on phase transitions in annular geometry, in both superfluids 
and superconductors. In such geometry, interplay of the symmetry breaking with the torus 

5 



topology of the system results in particularly dramatic predictions. Section 5 contains a 
an application of these considerations to the cosmological phase transitions. Conclusions 
are summed up in Section 6. 

2. SPONTANEOUS SYMMETRY BREAKING IN HIGH ENERGY AND 
IN LOW TEMPERATURE PHYSICS 

The ability to draw parallels between the high energy and condensed matter 
phenomena rests on the similarity of the behavior of free energy, especially in the vicinity 
of the phase transition. In particular, for the second-order phase transitions, potential 
contribution to the free energy is of the Landau- Ginzburg form: 

V{^) = a{T -T,)\^\'' +^-(3\^\\ (2.1) 

where is the order parameter or a field (see Fig. 1). The coefficients of the two forms 
have a well-prescribed dependence on the relative temperature: 

e = (T - Te)/Te . (2.2) 

Thus: 

a ~ a'e , (2.3) 

and 

/? = const. (2.4) 

In the field theories the form of the potential energy, Eq. (2.1), is often simply 
postulated (although it can be also justified by the Gaussian approximation at a finite 
temperature; Aitchison and Hey, 1982). In the low-temperature (or, more generally, 
condensed matter) context, it is, on the other hand, usually derived in the mean field 
approximation from the underlying microscopic theory of the system in question. For 
instance, in the case for superconductors, the so-called Gorkov equations provide a link 
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between the microscopic BCS theory and the Landau-Ginzburg theory (Gorkov, 1959). In 
any case, the parallels between the equilibrium aspects of the phase transitions in the field 
theories relevant to high-energy physics and in the effective field theories emerging in the 
mean field description of condensed matter system have been appreciated for quite some 
time (Aitchison and Hey, 1982, and references therein). 

Symmetry breaking arises when T < Tc, that is when the coefficient a in Eq. (2.1) 
becomes negative. Then the global minimum of the potential energy given by Eq. (2.1) 
moves from ip — (where it resides at T > T^) to the finite absolute value given by 



When (f is two dimensional instead of just two alternatives (+cr and — cr), degenerate 
minima form a circle of radius a. For 3-D </?, a sphere of radius a constitutes the degenerate 
manifold of the alternative true vacua. The depth of the minimum of a is given by 



The tension between the long range order which is supposed to set in below the critical 
temperature and the relatively short range over which the choice of the broken symmetry 
vacuum can be communicated and in a finite time determined by the rate at which the 
phase transition is taking place is responsible for creation of the topological defects. The 
tension between the long range and locality constitutes therefore the primary focus of the 
dynamics of nonequilibrium phase transitions which are the subject of this paper. 
Field Theories 

To consider consequences of the symmetry breaking in field theories, we supplement 
potential energy, Eq. (2.1), with a kinetic term and consider a Lagrangian L: 




(2.5) 



AV = V{0) - V{a) = a^/2p . 



(2.6) 



L(^) = d^^*d^'^ 



(2.7) 
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where <^ is a complex field. When a < 0, potential energy associated with Eq. (2.7) has 
a degenerate minimum — a circle at the radius a = ^/\a\Jp. Thus, the simplest (static, 
space- independent) solution of Eq. (2.7) is given by 

ip{x) = ae'^ , (2.8) 

where 6 is the fixed phase. Small perturbations around such a uniform solution can be 
considered. That is, one may investigate the behavior of 

(p{x, t) = (T+ {u{x, t) + iv{x, t)) /\/2 , (2.9) 

where u,v <^ a. Above, u is the perturbation of the absolute value of <f, and v is the 
perturbation of its phase. 

When L is expressed in terms of u and v defined by Eq. (2.9) and the constant terms 
are ignored. 



obtains. One can regard it as a Lagrangian for two coupled fields, u and v. Field u 
has a positive mass given by Pa^ = —a = and v — the so-called Goldstone boson — is 
massless. 

Let us now consider a more interesting static solution of Eq. (2.10). To do this, we 
begin by deriving a dimensionless version of the Lagrangian L, which obtains when cp, the 
distances, and the time are expressed in the "natural" units: 

<p^n = <p/a (2.11) 
f^g = f/^ (2.12) 
t^i=t/T (2.13) 

where a is given by Eq. (2.5), while ^ and r are the correlation length and relaxation time 
and are given by 

C = T = ^ (2.14) 



in the convenient cosmological unit system = c = A;^ = 1) we are employing in this 
section. Now a time-independent solution will have to satisfy 

V^=(|r7p-l)r7. (2.15) 

Apart from the trivial solution \r]\'^ = 1, Eq. (2.15) has an axisymmetric solution of the 
form: 

^/ = /(^')exp(m0) , (2.16) 

where (g, 0, z) are the cylindrical coordinates. Above, n must be a whole number — 
otherwise, r] could not be single valued. The radial part of the above solution is regular near 
the origin {f{g) — for ^ <^ 1) and approaches the equilibrium value at large distances 
(/^(^) ~ 1 — n^/ g"^ for ^ 1). The phase of this solution is 6 = ncj) on any g = const. > 
circle and is undefined for ^ = 0, where the amplitude of the broken symmetry phase 
decreases to zero. 

Equation (2.16) represents, of course, a string — topological defect — with the original 
symmetric vacuum locked out along the axis of symmetry by the broken symmetry phase 
in which a different vacuum is selected along the different radial directions in such a fashion 
that the transition to the minimum energy state, 77 = 1, is prohibited by an infinite energy 
barrier. The string tension — energy per unit length of the object defined by Eq. (2.16) — 
is logarithmically dependent on the cutoff in the distance. This solution was actually 
considered first by Ginzburg and Pitaevskii in the context of Landau- Ginzburg analysis 
of superfiuid helium (Ginzburg and Pitaevskii, 1958), and we shall soon return to its 
interpretation as a superfiuid vortex. 

An even more common and interesting field-theoretic model corresponds to a 
Lagrangian: 

L = —B^.B'^'' + [{d^ + ieA^) <f*] [{d^ - ieA^) <f] - a<f*^ - ^{<fi*<ff (2.17) 
with Afj, a massless gauge boson and 

= d^A, - d,A^ . (2.18) 
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While the Lagrangian given by Eq. (2.10) is invariant under a global gauge symmetry (i.e., 
redefining the phase 9 globally), Eq. (2.17) respects a more physically interesting local 
gauge transformation in which the phase is a function of space, that is 

^{x) e-*^(^V(a;) (2.19) 

and 

A^{x) ^ A^{x) - e-^d^e{x) . (2.20) 

When the expansion around the local minimum of |<^| = u we have implemented before 
for the case of the global gauge with the help of the substitution, Eq. (2.9), is carried out 
for the locally gauge invariant theory of Eq. (2.17), we will find 

^ z z (2.21) 

IBa'^u'^ - \/2e(jAf^dfj,v + . . . 

The surprising (and famous) result of symmetry breaking is the term involving A^A^ — it 
appears that the gauge vector field has acquired a mass (Aitchison and Hey, 1982). 

The excitations of the amplitude of (p around a are, of course, still massive. To see 
this, we can further simplify Eq. (2.21) by fixing the gauge so that 9{x) is the phase of the 
original complex field (p{x). Then Eq. (2.9) reduces to 

(p{x) ^ a + u{x)/V2 , (2.22) 

and the Lagrangian given by Eq. (2.21) becomes 



+ 



(2.23) 



In this form it is apparent that L describes interaction of the massive vector boson A'^ 
with a real scalar field m, the Higgs boson, with the mass given by 



ruH = P(T = \a\ . (2.24) 
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Vortex solution still exists (Nielsen and Olesen, 1973) for the Lagrangian given by 
Eq. (2.17), but its properties are now different than those given by the global string 
of Eq. (2.16). We shall consider it explicitly in a specific gauge below, in the course 
of the discussion of flux tubes in type II superconductors. Here, let us only anticipate 
one big difference: String tension associated with the local string is finite (rather than 
logarithmically divergent as was the case for global strings). 
Superfluid Helium 4 

The classic condensed matter example of a state with nontrivial broken symmetry 
phase is superfiuid He^. Phase diagram of helium is shown in Fig. 2. Superfiuid Hell 
occupies the low temperature (T < 2° K) low pressure {p < 25 atm) corner of the 
phase diagram. Astonishing properties of superfiuid helium (Tilley and Tilley, 1986) are 
associated with the Bose condensate of He^ atoms. Below the pressure dependent critical 
temperature (which is also often called T\ or "the A-line" because of the behavior of 
specific heat in its vicinity) , a temperature dependent fraction of atoms condenses out into 
a quantum fiuid described by a Bose condensate wave function 

* = |*|e'^ . (2.25) 

It is the phase ^ of ^' which is responsible for many of the intriguing properties of superfiuid 
helium. Some of them can be understood when we assume that obeys a Schrodinger 
equation of the form 

''^f = -|^^'*+''*- <2-2«> 

where m can be taken to be the mass of He^ atom, and /j, is the chemical potential — the 
energy gained by the system as a result of the addition of a single particle at constant 
volume and entropy. 

The equation for superfiuid fiow can be obtained from the Schrodinger equation (2.26) 
by computing the rate of change of where this probability density is now interpreted 
as the density of the superfiuid (Tilley and Tilley, 1986). Thus, for example, one can show 
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that the velocity v of the superfluid and the phase ^ of are connected with a simple 
equation: 

v=—Ve. (2.27) 

TO 

The formal connection between superfluid helium and our field-theoretic considerations of 
the previous section can be established when it is assumed that the free energy density 
can be expanded in powers of j^^p (which plays the role of the order parameter) and that 
it has a Landau- Ginzburg form: 

2 

(2.28) 

When the chemical potential /i is evaluated with the help of the above ansatz, the 
Schrodinger equation (2.26) we have written before acquires a form: 

zn— = -— V^* + a* + /3|*P* . (2.29) 
at 2m 

This relation is known as the Gross-Pitaevskii equation. It is generally treated with 
suspicion by the condensed matter theorists, who cite evidence that the so-called time- 
dependent Ginzburg-Landau (TDGL) equation (which has, on the left-hand side, a real 
rather than imaginary number multiplying the time derivative 9W / dt) is a better model for 
many systems (including superfluid He^) in the vicinity of the phase transition (Hohenberg 
and Halperin, 1977). 

Indeed, in the related case of superconductors TDGL equations have been rigorously 
derived from the microscopic theory (Gorkov, 1959). In that case the system is deflnitely 
overdamped, with the left-hand side of the equation for the order parameter dominated 
by r d^/dt where F is real. Although the microscopic theory of He^ superfluidity does 
not exist, the reason for overdamping (which in the superconductors is due to the coupling 
between the super and normal electrons) is certainly also present in superfluid helium. 
Thus, one can expect that approach to equilibrium will not be governed by Eq. (2.29), and 
this is indeed the case. The TDGL equation represents such relaxation quite accurately. 

12 



/(r) = a|*(f)|V||vE'(r)|V— V*(r) 



Nevertheless, Eq. (2.29) does have some interesting features. For example, the dispersion 
relation which obtains from the Gross-Pitaevskii equation does have both massless (a; ~ k) 
and massive (a; ~ A;^) components, although its form is not reminiscent (Bernstein and 
Dodelson, 1991) of the (second sound) phonon-roton dispersion relation which emerges 
only when the interatomic forces are taken into account (Vinen, 1969; Feynman, 1954). 

We shall not depend on the Gross-Pitaevskii equation (2.29) in analyzing details of 
the rapid phase transition (the subject of the next section). It is nevertheless good to write 
it down and use it to extract (with the help of the Landau-Ginzburg free energy) several 
important pieces of information which turn out to be largely independent of this (rather 
suspect) method of their derivation. The form of Eq. (2.29) differs from the equation of 
motion one would obtain from the globally gauge invariant Lagrangian, Eq. (2.7), only in 
that it has a first (rather than second) time derivative on the left-hand side. As a result, 
when Eq. (2.29) is re-scaled to assure a dimensionless form, 

ir,= -W^rj+{\v\^-l)rj, (2.30) 

the following three important quantities — relaxation time r, correlation length ^, and 
equilibrium density of the superfluid cr^ — are employed: 

r = n/\a\=To/\€\, (2.31) 

e = n/y/2m\a\ = Co/v^ , (2.32) 

and 

(7^ = -a/p . (2.33) 

Scalings with temperature follow when we adopt the usual ansatz [Eqs. (2.2)-(2.4)] for the 
parameters a and (3 (see Ginzburg and Pitaevskii, 1958, for numerical estimates of a' and 

/?): 

a = a'e, a' = 10"^^ [erg] ; (2.34) 

/3 = const. , /? ^ 4 • 10-^° [erg cm^] , (2.35) 
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where e = (T — Tg) /Tg, Eq. (2.2), is the relative temperature. 

When we restrict ourselves to the time-independent solutions of Eq. (2.30), we will 
be led back to Eq. (2.15) that we have already considered in the preceding subsection. 
In particular, the global string solution given by Eq. (2.16) for n = 1 is known as the 
Ginzburg-Pitaevskii vortex in the superfluid helium context (Ginzburg and Pitaevskii, 
1958) (see Fig. 3). The flow around the axis of symmetry is caused by the phase gradient, 
so that at a distance r away from the core of the vortex, its velocity is given by 

v = ---. (2.36) 
m r 

Vortices with the winding number n > 1 are unstable, as they are energetically more 

expensive than the n — 1 case. This can be seen by computing string tension which can 

be readily obtained from the kinetic energy of the flow for an arbitrary winding number n 

out to some cutoff L: ^ 

^ = / 27rrdr- (m|^'|^t;^) /2 

,2 . ... (2-37) 
n ■ • m 



Clearly, topological constraint (i.e., a certain winding number N within a closed path at 
some large radius) can be satisfied at a lesser energetic expense by N vortices with n = 1 
each rather than by a single vortex with n = N, for any N > 1. 

The string tension, Eq. (2.37), can be re- written in a more suggestive way: 

c; = [ne ■ {aV2f3)] •41n(^|) . (2.38) 

This equation [obtained from Eq. (2.37) in the case of n = 1] lends itself to an intuitively 
appealing interpretation: String tension is approximately proportional to the energy of the 
symmetric vacuum which occupies the core (of radius ^) of the vortex (and which is given 
by the term in square brackets). 

We also note the logarithmic dependence of <; on the large-scale cutoff L. [The small- 
scale cutoff at ^ appears naturally because of the absence of the superfluid condensate at 
small radii, where l^'p ~ , see discussion below Eq. (2.16).] 
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Existence of vortex lines was postulated nearly fifty years ago by Onsager (1949) and 
(independently, but somewhat later) by Feynmann (1954) to allow for superfluid rotation 
in spite of the requirement of single-valuedness of the Bose condensate wave-function. 
Vortex lines have since been studied both experimentally and theoretically (Tilley and 
Tilley, 1986; Donelly, 1991), although their possible origin in the "cosmological" scenario 
seems to have evaded attention of many of the low temperature scientists (see, for example, 
discussion of the vortex creation during the superfluid phase transition in Donelly, 1991, 
section 5.7). 

In addition to the vortex lines, we shall consider "cosmological" implications of 
persistent superfluid flows in an annular container. This phenomenon is again a 
consequence of the crucial role of the quantum phase of the Bose condensate in Hell: When 
the phase 9 increases by 27rn as one follows a closed path along the (large) circumference 
C = 27rr of the torus containing superfluid He'^, the velocity associated with the gradient 
of phase will be given by Eq. (2.27), 

, ?i , 27rn % n , , 

{-).-- = . 2.39 

m C m r 

Such flows are also stable (although for a somewhat different reason than the flow around 
the vortex line). Moreover, now n is essentially arbitrary [that is, until v reaches critical 
velocities at which superfluid flow begins to be dissipative (Tilley and Tilley, 1986)]. 

We end this subsection with a caveat: Landau- Ginzburg theory is only a qualitatively 
correct approximation for superfluid He^. Therefore, scaling relations for various quantities 
with the relative temperature e are somewhat different than one would infer from 
Eqs. (2.32)-(2.35). For instance, measurements show (Ahlers, 1976) that the correlation 
length ^ scales as , where = 2/3 and is consistent with the predictions of the 
renormalization group approach (sec, e.g., Goldenfeld, 1992). We shall pay attention to 
the consequences of such differences in what follows. Here, let us only note that the 
discrepancy between the Landau-Ginzburg approach and the behavior of superfluid He^ is 
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not unexpected. This is because, on the one hand, mean field theory which would allow 
one to apply the potential energy model of Eq. (2.21) to systems with atomic substructure 
is valid only when the correlation length ^ is much larger than the interatomic spacing a, 

^ > a . (2.40) 

On the other hand, considerations based on the Landau-Ginzburg model are valid only 
when the effects of thermal fluctuations can be neglected. This means that the inequality: 

Eg = ^^-AV> ^ksTc , (2.41) 

where AV = q;^/2/5, Eq. (2.6), known as the Ginzburg condition should be satisfied. 

In superfluid He^ the correlation length far away from T\ is of the order of 4-5 
Angstroms. Therefore Eq. (2.40) would demand e <^ 1. But (with the Landau-Ginzburg 
scaling ^ ~ l/->/|e[) the left-hand side of Eq. (2.41) scales as ■\/\e\^ and the inequality is 
satisfied only approximately 0.5°K below Tx. 

Hence, the two conditions, Eq. (2.40) and Eq. (2.41), cannot be simultaneously 
satisfied. As a result, the Landau-Ginzburg theory can be only applied as a qualitative 
guide to the behavior of superfluid helium (Tilley and Tilley, 1986). We shall return to 
the condition of Eq. (2.41) below, as it defines Ginzburg temperature: 

This relative temperature characterizes the demarcation line above which the potential 
difference between the minimum and the central peak is suflBciently small to make it easily 
traversable by the thermal fluctuations in the correlation length sized regions. 
S up er conduct or s 

Superconducting phase transition can be modeled by a Landau-Ginzburg mean fleld 
theoretic model with the free energy given by (Tilley and Tilley, 1986; Tinkham, 1985; 
Werthamer, 1969) 
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2 ^2 

+ — +Eo. (2.43) 

OTT 



Above, ^ is the order parameter — wave function of the Bose condensate of Cooper pairs 
with mass 2m and charge 2e — while A is the vector potential and 

S = V X 1 . (2.44) 

Energetic contributions due to the applied (external) field are incorporated in the constant 
Eq. 

Equation (2.43) can be justified rigorously by appealing to microphysics. The 
connection relies on the effective Gorkov equations (Gorkov, 1959) which can be derived 
from the BCS theory (Bardeen, Cooper, and Schrieffer, 1957) of superconductivity in the 
e <S 1 region. Thus, in contrast to superfluid He^, the Landau-Ginzburg theory is not just a 
phenomenological, qualitative approximation, but a microscopically justifiable quantitative 
tool (Tilley and Tilley, 1986; Tinkham, 1985). This fact will validate and simplify our 
discussion. On the other hand, instead of just one He^, the list of superconducting materials 
is rather long, even if we were to ignore the new high-T^ additions. Moreover, existence of 
the gauge field adds another characteristic length scale to the problem. 

Thus, by the same token as before, we can define the relaxation time scale, 

T = n/\a\^To/\e\ , (2.45) 
where tq can be derived from the Gorkov equation 

Correlation length: 



^ = h/^/^^\=^o/^M (2.48) 

characterizes variations of the order parameter while London's penetration depth A is 
proportional to the correlation length. 



A(T) = K ■ e(T) , 
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(2.49) 



where k is the temperature-independent constant which distinguishes between the two 
kinds of superconductors: Type I (where k < 1/V2) and type II (where k > 1/V2). 

For us this distinction is critically important, as the vortices can exist only in 
type II superconductors. There, the axially symmetric nontrivial solution which is a local 
minimum of Landau-Ginzburg free energy, Eq. (2.43), has the order parameter with the 
absolute value: 

hr 

(*/(7) ~tanhy , (2.50) 

where 6 is a constant close to unity. In the limit of extreme type II superconductors 
(k ^ 1), the induction varies as 

2^^' I A 



Bir) = ^,Ko{-) , (2.51) 



where Kq is a zero-order Hankel function of imaginary argument, and $o is a fliix quantum: 

$0 = hc/2e . (2.52) 

Qualitatively, B{r) falls off as e~^/^ at large radii (r > A) and can have significant values 
only in the shell ^ < r < A where it diverges as ln(A/r). This shell disappears in type I 
superconductors (for which ^ > v^A). Hence, type I superconductors do not allow for 
existence of vortices. 

String tension of vortex lines is finite, and it is given by 

<i^(^f\nK=^-A7r^^\nn, (2.53) 

where He is the so-called critical field (which, when applied externally, expels supercon- 
ducting phase from the sample). This formula is interesting because it demonstrates that 
the energetic cost of getting rid of the Bose condensate is given by the energy density 
associated with the symmetry breaking, AV of Eq. (2.6), so that 

AV = a^/2j3 = Hi /Stt . (2.54) 
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Thus, Eq. (2.53) implies that the string tension corresponds to the energy of the condensate 
displaced by the tube of radius r of normal metal times of proportionality factor of 4 • In 
Similarly, one can check that ^ is of the order of the energy of the field configuration 
associated with the flux line. 

As was the case for superfluids, one can also consider loops of superconductors in an 
empty space. Single- valuedness of the order parameter requires that the phase accumulated 
along the circumference of the loop must correspond to an integer multiple of 27r. Now, 
however — in contrast to the case of He^ — there are two sources of the phase gradient: In 
addition to the velocity v of the "superfluid" of Bose-condensed Cooper pairs, there is also 
the vector potential. Consequently, the quantization condition is 



In the limit in which the superconducting loop is wider than the Landau penetration depth 
A (so that there are closed paths C along which the first term containing v is negligible), 
Eq. (2.55) implies quantization of the flux in units of $o- This need not always be the 
case, as we shall find out while analyzing phase transitions in the annular geometry. Then 
A(r) becomes large, and sufficiently near Tc the resulting fiiix will be quantized in units 
smaller than ^o, with the velocity term making up the difference of the phase of the Bose 
condensate. 

The reliability of the Landau-Ginzburg theory for superconductors is due to the 
significant values of (typical ~ 1000 A, more than two orders of magnitude larger 
than in He^). Thus, the Ginzburg energy which is approximately given by 




(2.55) 
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(2.56) 



[see Eq. (2.54)] is typically much larger than the corresponding quantity for the same e 
in superfluid helium. To get an order of magnitude estimate of the relative Ginzburg 
temperature ec, we compute 



Thus, thermal fluctuations large enough to alter the configuration of the field are found 
only in the immediate vicinity of the critical temperature, even when a rather modest 
value for He (which can be ~ lO^Oe) and (which is typically closer to a few thousand 



It is perhaps useful to point out that the above estimates do not apply to the high 
temperature superconductors, which have small and much more modest critical fields 
He (Salamon, 1988). 
Other Systems 

The above two condensed matter phase transitions — the superfiuid He"^ and super- 
conductors — do not exhaust the list of the experimental possibilities, and the aim of this 
section is to point out — very briefly — that there are other systems in which topological 
defects in a nonconserved order parameter can be (or, indeed, have been) studied with the 
"cosmological" motivations in mind. 

Liquid crystals were already mentioned in this context. They were the first system 
in which a version of the experiment suggested for superfiuid He^ was carried out (Chuang 
et al, 1991; Bowick et al, 1994). Nematic liquid crystals consist of rod- like molecules (see 
for example de Gennes, 1974). Above the phase transition temperature these molecules 
are oriented randomly, but as the temperature falls they tend to align. It is not difficult 
to imagine how this local order may lead to topological defects. 

Formaly, the order parameter in the liquid crystal is given by the director field n(r), 
which is just the local orientation of the rods. By definition \n{r)\ = 1. Moreover, since 
the molecules can be fiipped by the angle tt without changing the physical configuration 




(2.57) 



Angstroms in type I superconductors) is adopted. 
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of the system, n{r) = —n{r). The vacuum manifold is therefore a two-sphere with the 
antipodal points identified. 

Topologically stable defects which are expected in such a symmetry breaking scheme 
include strings (also known as disclination lines) and a monopole. Rod - like molecules 
are arranged in a "hedgehog" configuration around a monopole. And there are several 
distinct kinds of linear defects. Rods can, for example, simply point away from the core 
of the string. In that case, when the core is circumnavigated the director field changes its 
orientation by 27r. But this is not the only possibility. Topologicaly stable defect can exist 
also when the director field changes its orientation by only tt, or by as much as Stt. This 
is because n{r) and —n{r) are identical. These topological defects can transform into each 
other in course of interactions. 

Phase transition into the nematic liquid crystal is of the first order. Following such a 
phase transition large densities of topological defects (especially the string-like disclination 
lines) were detected and studied (Chuang et al, 1991; Bowick et al, 1994). Disclination 
lines are visible under a microscope. Therefore, their scaling dynamics can be followed in 
detail. The disadvantage (from our point of view) is that the relevant phase transition is 
(weakly) first order. Thus, the most intriguing first stages of the development of the new 
broken symmetry phase can be understood in the relatively straightforward terms of the 
traditional nucleation paradigm, and do not shed any new light on the behavior of the 
order parameter in course of the second order critical dynamics. 

Superfluid He^ has even richer collection of topological defects. This is because the 
elementary particles which form the Bose condensate are pairs of He^ atoms, and, as such, 
have a nontrivial internal structure. This results in additional degrees of freedom which 
are accessible even at the very low (miliKelvin) temperatures in which He^ undergoes the 
superfluid phase transition. These additional degrees of freedom are associated with the 
non-zero spin of the Cooper pair (S=l) and with the orbital angular momentum of the 
Cooper pair around its center of mass. In the limit where the spin and angular momenta of 
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the Cooper pairs are strongly correlated, a liquid crystal like magnetic ordering will emerge. 
Thus, in addition to the broken gauge symmetry (present already in He^) rotational 
symmetry of both spin and orbital degrees of freedom can be broken. Consequently, 
the order parameter of superfluid He^ is described by macroscopic wave functions like in 
He^, but with nine complex amplitudes. Moreover, there are two different phases of He^ 
superfluid: The quasi-isotropic phase B, which corresponds to Cooper pairs with no total 
(that is, spin plus orbital) angular momentum, as well as the anisotropic phase A, for 
which total angular momentum is non-zero when projected along some direction. 

As a consequence of the rich structure of the order parameter, varied and interesting 
topological defects occur in both of the superfluid phases of He^ (Volovik, 1992). Moreover 
one He"^ version of the quench experiment has been considered already some time ago 
(Mineev, Salomaa, and Lounasmaa, 1987; Salomaa, 1987; Salomaa and Volovik, 1987), 
and, as we shall discuss shortly, another one has been recently carried out (Ruutu et al, 
1996). 

Helium 3 has both advantages and disadvantages when compared with He^. It is 
experimentaly more diflBcult to work with, and the variety of the topological defects may 
complicate the interpretation of the experimental results. Moreover, its correlation length 
is orders of magnitude larger than in He^, so the density of the defects is expected to be 
significantly less. On the other hand, large correlation length implies that the Landau- 
Ginzburg theory is a much better model, so the analogy with the field theories is better 
justified. Moreover, vortices can be detected "one by one" using nucler magnetic resonance 
along with some ither clever experimental techniques (Ruutu et ai, 1996). Hence, lower 
initial vortex line density may turn out not to be a problem. 

Last but not least, phase diagram of He^ allowes for a variety of quench trajectories, 
some of which raise the possibility of addressing other interesting issues such as the creation 
of topological defects in tuneably weak first order phase transformations, or even in a quick 
succession of the second and first order phase transitions. 
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3. QUENCHING OUT VORTEX LINES 

The analogy between the static, equihbrium properties of the field theories and 
condensed matter systems implies — as we have seen it in the preceding section — existence 
of similar kinds of topological defects. The purpose of this section is to employ this analogy 
to devise laboratory experiments which would allow one to study defect formation in the 
course of rapid phase transitions. I shall focus on the condensed matter systems with non- 
conserved order parameter only and consider defect formation starting with the simplest 
case — superfiuid helium. After a discussion of the freeze-out scenario, I shall go on to 
consider its consequences for the specific experiment already carried out by the Lancaster 
University group of Peter McClintock and his colleagues (Hendry et al, 1994; Hendry et al, 
1995). The theoretical analysis will be repeated for the even more interesting (local gauge) 
case of a rapid phase transition to a superconductor, although it will be pointed out that 
the main obstacle for a successful study of a superconducting analog of the cosmological 
scenario is the difficulty of rapid detection of the vortices (fiux lines). Thus, discussion of 
a more experimentally attractive (although more theoretically challenging) version of the 
superconducting experiment — which involves locking out magnetic flux in a loop — will be 
delayed to Section 4. 
Quench Into Bulk Superfluid 

Superfluid He^ has a number of experimental advantages which make it — in spite of 
the theoretical complications we have hinted at in the previous section — perhaps an ideal 
laboratory to study dynamics of phase transitions in systems with global symmetry. Thus, 
for example, one can obtain samples of He"^ of almost unparalleled purity. Furthermore, 
He'^ can accommodate only one kind of defect — vortex lines — simplifying the analysis. 
Moreover, detection of vortex lines — while not as straightforward as it is for defects in 
liquid crystals, where they are visible almost with a naked eye (Chuang et a/., 1991; 
Bowick et ai, 1994; Yurke, 1995) — is also relatively simple, as it relies on the well-tested 
techniques of second sound attenuation (Donnelly, 1991). We have already mentioned 
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that the superfluid transition is of the second order. And last, but definitely not least, 
phase transition in He^ can be induced on a dynamical time scale by rapid decompression. 
Therefore, the critical temperature Tx will be reached throughout the bulk of the He^ 
sample on a very short time scale limited by the first sound travel time (~ 220 m/s in 
the vicinity of the critical temperature T\). This is orders of magnitude larger than the 
(second sound) velocity defined by the ratio of the correlation length and relaxation time 
scale: 



which limits the speed with which different regions of the emerging superfiuid become 
correlated. Above, io/rQ evaluates to approximately 70 m/s but second sound reaches this 
velocity only quite far from Tx. In the vicinity of Tx s is much less, so that the "second 
sound horizon" can be quite a bit smaller than the size of the sample. Second sound gives 
the rate of propagation of perturbations of the density perturbations of the superfluid. 
Hence, it also limits the rate at which the order parameter will be able to adjust. In this 
sense, one is able to reproduce the "acausal" nature of the cosmological phase transitions in 
the superfluid He"^ better than in the other systems, especially where the transport of heat 
is involved (as is the case in liquid crystals and will likely be the case in superconductors). 

The schematic quench trajectory (which proceeds along the isentrope) is shown in 
Fig. 4. Figure 5 illustrates the transformation of the effective (Landau-Ginzburg) free 
energy, as well as its hypothetical effect on the conflgurations of the order parameter. The 
aim of this section is to estimate the density of the vortex lines left behind by the quench. 
To do this, we shall rely on the critical slowing down — the behavior of the relaxation time 
scale T in the vicinity of Tx- As we have seen in the previous section (and as is conflrmed by 
the experiments), r ~ l/|e|. Thus, in the course of the quench the time scale on which the 
order parameter can adjust to the new thermodynamic parameters (especially to the new 
value of e!) is becoming very long in the vicinity of the critical temperature. As a result, 
two regimes can be distinguished: (i) Sufficiently far from Tx the relaxation time scale r 




(3.1) 
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is much smaller than the time on which the quench is proceeding. In this adiabatic regime 
order parameter will be characterized by an equilibrium configuration with the correlation 
length ^ determined (to an excellent approximation) by the instantaneous value of e. By 
contrast, very near Tx the equilibrium relaxation time scale will be much larger than the 
time spent by the system with the corresponding value of e. As a consequence, we can 
define: (ii) The impulse region, in which r is so large (because of the critical slowing down) 
that the configuration of the order parameter will be in efi^ect immobilized on the timescale 
of interest. 

The boundary between these two regimes will occur at the freeze-out time t. We 
can compute it by assuming that, in the vicinity of Tx, the relative temperature e is 
approximately proportional to time: 

e = t/TQ . (3.2) 

Quench timescale tq can be controlled by the rate at which the pressure is lowered. Critical 
temperature is reached at t = 0. 

Freeze-out time is then set by the equality: 

T{i) = i , (3.3) 

which, given Eq. (3.2) above, as well as r = To/|e|, yields: 

i=,/f^. (3.4) 

Consequently, the transition between the adiabatic and impulse regimes occurs (twice) 
during the quench at the relative temperature: 

e = e{i) = ^tq/tq . (3.5) 

Here, tq ~ ^/ck' ~ 10~^^ s, with more careful estimates yielding tq = 8.5 • 10~^^ s. The 
fluctuating configuration of the order parameter is frozen out at e > 0, above Tx- When the 
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dynamics are "restarted" again with e < 0, the parameter a is also negative, which means 
that the symmetry breaking has aheady occurred and topological defects are "frozen out." 

The characteristic correlation length ^ corresponding to t and e above will decide the 
(initial) density of the vortex lines through the approximate formula: 

4 = i~r', (3.6) 

which is based on the idea that typically a piece of the vortex line will fall within a ^-sized 
domain (Kibble, 1976; 1980; Vilenkin, 1981). 

Correlation length diverges near the phase transition temperatures: 

^(e) = wkr- (3.7) 

In the mean field theory v = \/2 [see Eq. (2.32)] and = 5.6 A provide an acceptable fit to 
the data. However, experiments seem to point to a better accord with the renormalization 
group prediction for v — 2/2> and a somewhat different = 4 A to go with it (Ahlers, 1976; 
Goldenfeld, 1992). Consequently, the estimated value of the initial vortex line density is 
given by 

4-eo-^(^)^ (3-8) 

The first obvious remark is that — because of the smallness of — even for relatively slow 
quench time scales the anticipated density of vortex lines is enormous. The second comment 
is that this very large initial vorticity is likely to disappear on a rather rapid time scale. 
This is because the density of vortex lines decays approximately as (Vinen, 1957): 

I = -X-^' - -li' , (3.9) 
at m 

where the Vinen parameter x is a dimensionless constant. Thus, I will decrease rather 
quickly, with the time-dependent density given by 

£ = ^° . (3.10) 
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This process is also reminiscent of the dynamics of cosmic string network, although absence 
of cosmological expansion, the differences in the equation of motion of vortex lines, etc., 
make this a rather distant analogy. 

Our derivation of the preferred scale ^ and the resulting initial density of topological 
defects can be briefly summed up by noting that we are looking for the instant when the 
quench becomes effectively instantaneous. The correlation length at that instant will be 
obviously "frozen out" and will set the initial scale of the vortex line network in the regime 
where the symmetry of the order parameter is broken. 

There are two other ways of deriving the estimate of the characteristic scale in the 
course of second-order phase transformations which appeal to somewhat different physical 
input but lead to the same estimate of ^. We shall consider them now as they may be of 
heuristic value to the readers. 

To describe the first of them, we compare the velocity with which the correlation 
length would have to increase to always maintain its equilibrium value, 

v^ = ^ = id^/de)e , (3.11) 

with the velocity of the second sound or, more generally, with the speed of propagation of 
the perturbations of the order parameter, 

s = e(6)/r(6) . (3.12) 

Simple algebra shows that would have to increase faster than s when e is smaller than 
the characteristic value e which satisfies equation: 

v^{e) = s{e), (3.13) 



or 

eoe-^/V2Tg ^ ^oe'/^ro ■ (3.14) 

This last equality leads to Eq. (3.5), providing we ignore numerical corrections of order 
unity. The physical implication of this way of obtaining e and ^ are straightforward: ^ 
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is as large as the dynamics of the order parameter (characterized by s) and the quench 
timescale tq allow it to become. 

Freezeout temperature e can be also obtained by comparing ^ with the size of the sonic 
horizon h — the distance over which perturbations of the order parameter can propagate in 
the course of a quench. This distance is set by 





The system will emerge in the broken symmetry phase with the correlation length no larger 
than the sonic horizon. This will be associated with the relative temperature given by 



which, again, after some straightforward algebra, leads (approximately) to Eq. (3.5). 

The three above arguments lead to the same conclusion. This is because they are 
essentially equivalent. The only difference between them is the manner in which they 
all appeal to the same concept of two different regimes of the quench — (i) slow, nearly 
adiabatic for |e| > e and (ii) almost instantaneous, nearly impulse for |e| < e — by 
respectively comparing timescales, velocities, and distances characterizing dynamics of 
the order parameter in the vicinity of the critical temperature Tc- 

It may perhaps also be useful to note that the Vinen equation, Eq. (3.9), can be 
reexpressed as an equation for the coherence scale of the order parameter phase, D, by 
simply defining 



D{t) is the size of the domains over which the order parameter has similar phase. Using 
Eqs. (3.9) and (3.17), one easily obtains 



t 




(3.15) 



(3.16) 



D{t)^l/^. 



(3.17) 



dD 
2— 



dt 



l/D. 



(3.18) 
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This is, of course, readily transformed into 



d 
di 



(3.19) 



which immediately yields the solution 



D{t) = D{0) ■ . 



(3.20) 



This rewriting of the equation for the decay of vortex lines in superfluid helium in the 
language of domain size growth offers a different perspective on the evolution of the order 
parameter in the post-quench period. It also establishes a connection between the evolution 
of the vortex line tangle in He^ (see Donnelly, 1991) and the more general problem of 
evolution of the domain size in the systems with nonconscrved order parameter. There, on 
the basis of quite general considerations, one can establish Eq. (3.20) [except for the case 
of vortices/monopoles in two dimensions, when D{t) ~ (t/lnt)^^^] (Bray, 1994). 
Activation Mechanism 

The process we have ignored so far in the discussion — thermal activation — was initially 
believed to be essential in generating topological defects (Kibble, 1976). As we shall see 
below, this expectation does not seem to be borne out in the case of the superfluid phase 
transition in He^ [nor was it anticipated in the original experimental proposal (Zurek, 
1984, 1985, 1993)]. To see why vortex line density is unlikely to reach values mandated 
by the thermal activation processes at the Ginzburg temperature, we consider the initial 
vortex line network laid down at t by the freezeout of the order parameter. Approximately 
70% of the vortex lines will be in a form of a long string, with only the remaining 30% 
in the form of loops, the smallest of which have size ^. This follows from the numerical 
simulations (Vachaspati and Vilenkin, 1984), and approximately coincides with the fraction 
of nonintersecting random walks, although the precise numbers depend on the lattice as 
well as on the distribution of the sizes of domains in a way which is not yet understood 
(Hindmarsh and Kibble, 1995). Let us now suppose that — as will be typically the case — 
we have stopped the quench at some temperature T below the freeze-out temperature 
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T = so that the original scale ^ of the defect network substantialy exceeds ^{T). 
Thermal fluctuations will generate loops of size ^, which will be typically smaller (much 
smaller) than the pre-existing structures. In particular, the long string can be modified 
due to the small-scale warping by the addition of perturbations on the scale ^; but this 
process has an essentially diffusive character, so it will proceed rather slowly. By contrast, 
the long string has a total energy which is much larger than the Ginzburg energy Eq, 
Eq. (2.41): 

Ev = Eg- ■ V « {^/i) ' • («V2/3) • V . (3.21) 

Above, V is the volume of the container. This equation is easy to understand — it can be 
readily rewritten as 

= £ • • (aV2/?) ~ £ • ? , (3.22) 

where £, stands for the total length of the string. Clearly, for any macroscopic-sized volume 
V of the superfiuid, energy Ey will be much in excess of the thermal energy. Consequently, 
while thermal activation energy may suffice to generate small-scale ringlets of the vortex 
line, it is unlikely to elongate the string. Indeed, if anything, the long string may begin 
to straighten before the quench is completed. This process will counteract — and will most 
likely even outweigh — thermal activation events which would act to increase vortex line 
density. 

Consequently, it is plausible that the dominant vortex line formation process in the 
course of a quench will be the freezeout of the order parameter (Zurek, 1984, 1985, 1993). 

Comparison with the He^ Experiment 

In the recent experimental implementation (Hendry et al., 1994, Hendry et al., 1995) 
of the cosmologically inspired vortex line formation scenario (Zurek, 1984, 1985, 1993), 
McClintock and his colleagues of the Lancaster University have carried out the quench 
through the A-line, approximately along the isentropic trajectories shown in Fig. 6. As 
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anticipated in the "freezout" scenario of defect formation, this led to a copious defect 
production with the initial line density, estimated to be 

li > 10^^ m-2 = 10^ cm-2 (3.23) 

in the quench which crossed the A-line (Hendry et al., 1995). This is certainly a very 
large vortex line density. If a comparable density were to be generated by rotation of the 
(~ 1 cm diameter) sample they have used, angular velocities of ~ 4-10^ radians per second 
would be required. 

It is possible to estimate the quench time scale tq from the change of relative 
temperature (Ae ~ 0.1) and from the time (At ~ 3 ms) it took to complete the pressure 
drop: 

TQ = At/Ae = 30 ms . (3.24) 

Using this estimate of the quench time scale in the equation (3.8) we have just derived, we 
are led to predict 

iLoim-''] ~ 3 • lO^V (t-q/100 ms)'/^ , (3.25) 
when Landau- Ginzburg scaling is employed. A somewhat more modest 

£flG[m-2] ~ 1.2 ■ lO^V (rg/lOOms)'/^ (3.26) 

obtains for the renormalization group scaling (Zurek, 1995). As it was already pointed 
out in the commentary (Zurek, 1994) on Hendry et al. (1994) [where the estimate was 
evaluated from Eq. (3.26) for tq = 30 ms and compared with the then available lower 
bound of £j > 10^^ m-^], the agreement between the (order-of-magnitude) estimates offered 
by the theory and dramatic experimental results is impressive. Of course, more could be 
done on both theoretical and experimental fronts. Indeed, progress on the experimental 
front may be more rapid (for example, if the proposals listed in Hendry et al., 1995, are 
implemented) than in the theoretical treatment of the problem (since microphysical theory 
of superfluidity in He^ is still missing) . 
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A few remarks on both the experimental side (where I will largely echo discussion of 
Hendry et al., 1995) as well as on theoretical aspects of the treatment are nevertheless in 
order. In addition to: (1) the quench through the A- line, McClintock and his colleagues 
have carried out two other types of quenches; (2) quenches far away from T\; and 
(3) quenches starting just below (few milli-Kelvins) the phase transition temperature. 
As we have already noted, copious production of vortex lines was detected in type (1) 
quenches. By contrast, no detectable vortex line production was detected in quenches 
of type (2) which never approached T\. This is reassuring, as it shows that vortex line 
creation can occur only where the transition crosses the A-line and freezes-out the pre- 
existing fluctuations of the order parameter. However, quenches of type (3) did create 
detectable vorticity, although in the amounts which are significantly lower than for the 
A-line crossing quenches of type (1). 

This is an intriguing observation. At first, one might be tempted to appeal to thermally 
generated seed vortices which could exist in the "Ginzburg" regime, but — for the reasons I 
have already outlined — it seems unlikely that long (i.e., more than a few correlation length 
in diameter) vortex lines may exist a few milli- Kelvin below T\. 

If vortices detected following the type (3) quench from just below the A-line are not 
generated by thermal fluctuations alone, what else can they be? And why are they detected 
in type (3) quenches but do not appear when the quench starts far from the A-line? I believe 
the crucial clue to the interpretation of this phenomenon may be the remark of Hendry et al. 
(1995), which attributes appearance of these extraneous vortices to the flows generated in 
He^ by the pressure quench. The difference in the generated vorticity in type (2) and in 
type (3) quenches may be due to the different cost — measured in the vortex line tension 
Eqs. (2.37) and (2.38) — of creating topological defects. Thus, if the same energy AE'stir 
were available in the stirred-up superfluid. 



^STiR ~ AEstir/? 
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(3.27) 



of vortex line length could be created. This simple reasoning would predict that the 
stirred-up vortex line density should scale (to the leading order) as 

^STiR - , (3.28) 

where we have recognized that the vortex line string tension q falls off with pg (and where 
we have ignored logarithmic terms; if they were taken into account, one would be led to 
the equation: ^stir • Ps ~ (~ Ih^stir^^) ^, where both the superfluid density and the 
correlation length ^ are evaluated near the initial point on the quench trajectory). 

This simple discussion ignores many effects which could be relevant, but it does have 
an advantage of being testable: Quenches along the trajectories which differ only in the 
initial, near-Tx points, should allow one to verify (or prove wrong) Eq. (3.28) above, 
perhaps with the logarithmic corrections. 

This discussion leads us to one further remark about the vortex line density generated 
in quenches which do cross A-line. It is conceivable that the vorticity generated by the 
cosmological mechanism in the course of the quench may be further amplified by the 
inadvertent stirring. If this were indeed the case [and discussion in Hendry et al. (1995) 
does admit such possibility] , then the direct comparison of the predictions of our discussion 
[e.g., Eq. (3.8)] and of the experimental results must be taken with a grain of salt. 

It would be therefore quite important to have — in addition to the lower bound given 
by Hendry et al., 1995 [Eq. (3.23)] — to also have an experimental upper bound. In view 
of the rapid evolution of the vortex line density [Eq. (3.9)], it might be useful to attempt 
setting a thermodynamic upper limit on After all, decay of the vortex line density will 
convert string tension energy into the heat deposited in the superfluid, which should allow 
one — especially for li ^ 10^^ m~^ — to set limits on the initial £i, simply by monitoring 
the temperature of the superfluid. The anticipated total temperature increase as a result 
of decay of the vortex line density I can be estimated: 

AT = Iq/Q . (3.29) 
33 



Above, Q is the appropriate specific heat of He^ at the temperature at which the decay 
of the network is taking place. Whether detecting corresponding temperature differences 
resulting from Eq. (3.19) may be possible on the scale rapid enough and with the accuracy 
sufficient to make the results useful is, of course, a rather difficult experimental question.* 

If one were certain that the only way in which the quench could deposit heat in the 
superfluid was by dissipation of the string tension energy locked in the vortex line network 
generated in the course of the quench, one could amplify the effect by traversing the A-line 
many times. Then the accumulated heat would be due to the small departures from the 
isentrope which are caused by depositing some of the energy in the ordered configurations 
of the Bose condensate, which — following each traversal of the A-line — are dissipated in 
the destruction of the vortex line network. Among the assumptions which are made here 
are an excellent isolation of the sample as well as a sensible guess that the energy locked in 
the vortex network after each quench comes from the mechanical energy available in this 
nonequilibrium process rather than from thermal fluctuations (which, again, relates to the 
difference between freeze-out and thermal activation). 

In summary, one should stress that all of the caveats listed above do not weaken the 
main conclusion of Hendry et al. (1994 and 1995). Very abundant vortex line creation 

* Using Eq. (3.29) above and a remark in Section 3 of Hendry et al. (1995) that "Fast 
expansions seem to follow the quasistatic ones to good approximation, but, from a common 
starting temperature, reach a final temperature which is typically several mK higher" we 
can set an upper limit on t. 

h5[m-^] < 0.7AT[mK] . 

Thus £i5 ^1 (or I ^lO"*^^ m~^) is a plausible upper limit on the vortex line density. 
Moreover in veiw of the inadvertant "stirring" in the course of the quench reported by 
Hendry et al. (1995), it might be a generous overestimate of the vortex line density 
attributable to "pure" quench. 
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is observed in quenches which cross the A-hne. Estimates based on the apphcation of 
the cosmological scenario of topological defect creation are consistent with the detected 
vortex line density, providing that one relies on the idea of the freeze-out time t rather 
than on thermal activation. Indeed, the final points of all the quenches are well within the 
Ginzburg region of He^, which would imply — if one were to rely on the vortex line creation 
by thermal activation — that the final vortex line density should not decay but, rather, 
that it should reach values set by the correlation length associated with the end-of-quench 
point. This is clearly not the case. Further experimental study of the freeze-out mechanism 
should involve limiting the contribution of "stirring up" of the superfiuid in the course of 
the quench as well as varying the quench time scale to verify Eq. (3.8). This course of 
research has been already mapped out by the Lancaster University group (Hendry et al, 
1995). 

Quenching out Defects in Superfiuid He^ 

The first demonstration of defects formation in course of the phase transition into the B 
phase of superfiuid He^ has been just reported by the Helsinki group and their collaborators 
(Ruutu et al., 1996). The transition ocurrs in a cigar-shaped, 100 //m sized volume in a 
rotating container of the B- phase of He^, which is heated above the thermal temperature 
by the decay products of a thermal neutron. As this volume cools and re-enters the broken 
symmetry phase, string-like defects form with the density which is consistent with the 
freezeout distance ^ = Co(t'q/to)^/^ of approximately 1 //m. Ginzburg relative temperature 
for the superfiuid He^ is very small, and corresponds to the size of the coherence length 
much larger than the volume of the "cigar" in which the phase transition occurs. Thus, 
again, the result of the experiment is inconsistent with the activation mechanism, but 
appears to conform with the freezeout scenario. 

The superfiuid is slowly rotating. As a result, sufiiciently large vortex loops generated 
in the small volume that underwent the phase transition are stretched by the Magnus force. 
This eventualy tends to transform them into rectilinear vortex lines parallel to the axis of 
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rotation. They are subsequently pulled into the center of the container, where they can 
be detected by NMR. The distribution of loop sizes obtained in this manner is consistent 
with the scale invariant distribution (see, e. g., Vachaspati and Vilenkin, 1984). 

This very exciting experimental development complements work on superfluid He^ 
(where the density of vortex lines is typically much higher, but where they are also harder 
to detect). It also opens up the possibility of pressure quenches into either of the He^ 
phases. In particular, phase diagram of He^ allowes for a variety of quench trajectories 
(see Fig. 7), some of which raise the possibility of addressing other interesting issues such 
as the creation of topological defects in tuneably weak first order phase transformations, 
or even in a quick succession of the second and first order phase transitions. 
Quench Into Bulk Superconductors 

Superconductors of type II have a Landau-Ginzburg form of the free energy, Eq. (2.43), 
and a negative surface energy (so that they can accommodate vortex lines). We can 
therefore repeat the argument we have put forward for the superfluid phase transition to 
compute the density of vortex (flux) lines which obtains as a result of a rapid quench. 
Again, we shall rely on the freeze-out time scale t and the corresponding correlation 
length of the order parameter ^ = ^{t). The differences between the case of He^ 
and superconductors will be not in the key ideas behind the estimates of the vortex 
line density but, rather, in the different experimental requirements, which will make 
superconducting bulk quenches more difficult to carry out. Moreover, the consequences of 
the superconducting quenches may not be easy to investigate (although, as we shall see in 
the next section, loop geometry may help!). 

As before, we imagine the quench proceeding on a time scale tq, with the distance 
away from the phase transition (measured by the relative temperature e) changing as 
e{t) = t/rq, Eq. (3.2). As the relative temperature varies, the relaxation time scale is also 
changing, r = To/|e|, so that sufficiently near the phase transition instant (where e = 0), 
dynamics of the order parameter becomes too sluggish to adjust the correlation length ^ 
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to the equilibrium value given by Co\/R' (2-48). The switch from the adiabatic to the 
impulse regime will happen (as it was the case for superfluids) at the freeze-out instant: 

i = V^oTQ , (3.30) 

where tq is the characteristic time given by Eq. (2.47). Consequently, we can evaluate 

£^ 1.225^tq/TcH , (3.31) 

where tq is in seconds and the critical temperature is expressed in Kelvins. The frozen-out 
correlation length is then given by 

i = 10-2 (Co/lOOO A) T^/"[cm] (3.32) 

and corresponds to the critical relative temperature: 

e = e(i) = 1.225 • IQ-^/y/T^ . (3.33) 

This e is rather small: Therefore, it might not be easy to go through the phase transition 
both sufficiently quickly and sufficiently uniformly to validate the standard prediction 

£ = 1/|2 = 10^(1000 A/eo)Vg/'[cm-2] . (3.34) 

The anticipated initial Aux line density is therefore much smaller than was the case for 
superfluids. The decrease in estimated £ is mainly due to the increase of typical ^o, which, 
for a typical superconductor, is two orders of magnitude bigger than for the superfluid 
He^. One could, of course, increase anticipated initial £ by selecting to work with high- 
temperature superconductors (which have much smaller ~ 10 A). This strategy may 
be worth pursuing, although the switch to high-Tc materials is likely to be accompanied 
by a specific set of problems (but, possibly, also by some advantages). In particular, the 
Landau-Ginzburg theory may not be as accurate for the reasons we have explored before. 
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Moreover, our estimates of tq will not apply, as they are based on the microscopic theory 
(BCS) which does not apply to high-Tg superconductors. 

Let us now return to the experimental complications. To begin with, the pressure 
quench will not result in as big a change of critical temperature as is the case for T\ in 
He"*. There is, of course, some sensitivity of critical temperature in superconductors on 
applied pressure (see, e.g., Lynton, 1963), but the effect is rather small and may be difficult 
to exploit experimentally. Critical temperature in high-Tg material is more sensitive to 
pressure (Neumeier and Zimmerman, 1993; Neumeier, 1994). This may be another reason 
(in addition to smaller ^o) which could make them interesting for our purposes. 

Temperature quench — rapid cooling of the sample — may be nevertheless the most 
practical way to precipitate the phase transition. Unfortunately, this is a diffusive process. 
Therefore, it is likely to be relatively slow. Moreover, heat transport is (of course) driven 
by temperature gradients. Thus, limiting the gradient of the relative temperature to much 
less than the value set by the ratio: 

9 = e/i , (3.35) 
may not be an easy task. Gradient g evaluates to: 

g ^ 1.225 • 10-''(1000 A/^o)7;"^/Vg^/'^[cm-i] , (3.36) 

when Gorkov equations (Gorkov, 1959; Tinkham, 1985) apply. Again, one may be tempted 
to consider high-Tg materials, but perhaps a more practical suggestion would be to work 
with thin layers of superconductors. This 2-D strategy allows one to transport the heat 
in the direction perpendicular to the plane in which the superconducting state will form. 
Thus, one may be able to violate the constraint imposed by Eqs. (3.35)-(3.36) without 
invalidating the reasoning which has led to the estimate of ^, Eq. (3.32). This two- 
dimensional sample strategy may have one additional advantage: The interior of the bulk 
superconductor is difficult to probe. However, flux lines can be detected when they emerge 
from the sample. In a sample with a thickness of the order of ^ (~ 10~^ cm) or even 
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somewhat larger, all of the vortex lines can be "counted" where they pierce the surface. 
Moreover, using 2-D geometry may help in slowing down annihilation of vortex lines as 
one may be able to supply in a controlled fashion sufficiently many pinning sites to prevent 
or at least significantly impede vortex line migration. 

The estimate of the vortex line density is now — in the 2-D case — given of course by 
Eq. (3.24) but with a somewhat revised interpretation. We are led to expect a flux line to 
emerge from within the area of the order of ^. Thus, the surface density of the vortex line 
endings is given by 

(3.36) 

Scattering of neutrons from a sample undergoing a superconducting phase transition 
may be a way to peer inside a truly 3-D fragment of material. Neutrons have a magnetic 
moment and will scatter from the network of flux-filled vortex lines, which should in 
principle allow one to deduce vortex line density and perhaps even get some information 
about their geometry. 

It has also been suggested (Rudaz, Srivastava, and Varma, 1994) that an electron 
microscope could be used to image individual flux lines, although it is harder to imagine 
how one could use it to track a network of lines in the course of the phase transition. 

So far we have ignored one more obvious complication. Superconducting vortices 
are associated with the gauge field (B) which has to be screened in the experiments. 
Moreover, the relative importance of the order parameter and of the gauge field in the 
phase transitions with the local gauge symmetry breaking is not entirely obvious (Rudaz 
and Srivastava, 1993; Kibble and Vilenkin, 1995). This last remark may be regarded more 
as a motivating factor (rather than as a reason for discouragement): There are real, deep, 
and interesting questions and scenarios of defect formation in theories with local gauge 
symmetries. It is however also conceivable that transport of the field through the material 
on the verge of becoming a superconductor may complicate matters in a more mundane 
manner. 
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4. COSMOLOGICAL EXPERIMENTS IN ANNULAR GEOMETRY 

Let us consider a system undergoing a symmetry-breaking phase transition. In the 
vicinity of the freeze-out instant, t, the fate of its order parameter is sealed: Phase 
differences in the neighboring domains are set. Hence, the distribution of the topological 
defects is also essentially decided. Each vertex defined by the domains with the phase 
differences which add to ~ 27r will become a section of a topological defect — a superfluid 
vortex, a flux line in a type II superconductor, a disclination line in a liquid crystal, or a 
cosmic string. This phase difference can be computed by adding the phases along a path 
which circumnavigates the vortex and which has a radius comparable to the size of the 
frozen-out domains. 

The same postulate of the independent choice of the phase in independent domains 
also allows one to compute the total phase difference along a path with much larger 
circumference C, which traverses many (jV) domains. The total phase difference along 
such a path is given by 

Mc^VM=\[c/1. (4.1) 

Already in the bulk experiments we have discussed so far this statement is of some interest, 
as it implies that the vortex lines will tend to be more anticorrelated than would be the 
case if their directions were assigned at random. This is easy to see. For, when the vortex 
line orientations are selected randomly, the net circulation along a path defined by a circle 
of size r (area 7rr^) would increase with the square root of the number of strings enclosed, 
and hence 

~ ^Jirr^/e ■ (4.2) 

This implies increasing with the circumference rather than with the square root of the 
circumference [Eq. (4.1)], as the more careful analysis implies. This discrepancy between 
the two predictions, Eqs. (4.1) and (4.2), can be, of course, readily settled in favor of 
Eq. (4.1), not just by the discussion we have already carried out but by an additional 
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observation (Vachaspati and Vilenkin, 1984) that the string network generated in the 
course of a rapid phase transition will consist (~ 70%) of a long string which will have to 
cross any surface spanned by C in opposite directions more often than a random collection 
of vectors. The rest (30%) of the length is due to (mostly small) loops which will typically 
not contribute to A^c unless C passes through their centers (i.e., unless they are "strung" 
on C like beads). 

The purpose of this section is to discuss how Eq. (4.1) can be tested directly in both 
superfluids and in superconductors. 
Superfluid in an Annulus 

Consider quench in an annular container of superfluid He^ (see Fig 8). The phase 
difference locked out in the course of the quench is given by Eq. (4.1) with 

e = wi^r, (4.3) 

where e is given by ^toTq, Eq. (3.5), and v = 2/3 in the renormalization group treatment 
of He^. Therefore, in accord with the equation (2.27) which relates the phase gradient with 
the velocity of the superfluid, one expects (after the quench) a flow in a random direction 

with the velocity of 

. = A.i,/S = A/X. (4.4) 

This is obviously an intriguing prediction, especially since it evaluates to measurable 
velocities. To estimate the resulting v, we can further calculate 

Above, To ~ 8.5 • 10~^^ s, and C is measured in centimeters. Thus, for tq of the order of 
milliseconds and u = 2/3 one obtains 

V ~ 0.4 (tqH)"^^^ /VC[cm] [mm/s] . (4.6) 

This is certainly a macroscopic (if not very large) velocity. Moreover, Eq. (4.6) is likely 
to be an underestimate, since Eq. (4.1) contains an implicit assumption that the phase 
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difference between the domains separated by ^ is approximately a radian, while the more 
usual assumption would be to set that difference at n or 27r/3, which would increase the 
estimate given by Eq. (4.6) above 1 mm/s. A still higher estimate of the quench-generated 
velocity (probably of the order of cm/s) would follow if we used a lower bound on the 
quench-generated vorticity implied by the experimental results of the Lancaster group 
(Hendry et al, 1995) to infer the relevant ^. 

In any case, locking out a finite "phase around the loop" as a result of a rapid 
symmetry-breaking phase transformation is the essence of the cosmological scenario for 
the generation of topological defects. Performing the experiment in an annulus is bound 
to be a challenge, but the dramatic prediction we have outlined above makes it a worthwhile 
effort. Moreover, the locked-out phase in the loop has one additional advantage over the 
vortex creation in the bulk. The resulting winding number 71^: 

riu, = A^c/27r (4.7) 

is [to an excellent approximation (Tilley and Tilley, 1986)] independent of time in the 
superfluid phase. Thus, quench can set up a persistent superflow with a macroscopic, 
measurable velocity. This velocity is directly related to the frozen-out domain size ^ 
through Eq. (4.4) and can be — in contrast to the rapidly decaying vortex line density £ 
[see Eqs. (3.8)-(3.10)] — measured at leisure long after the quench. 

This discussion brings up a few questions. Let us start with an apparent paradox: 
How can a system which did not rotate before the phase transition acquire a finite 
angular momentum [corresponding to the velocity of f , Eqs. (4.4)-(4.6)] as a result of 
an axisymmetric pressure quench? A facile answer to this question would be to simply 
point out that the considerations concerning the order parameter of the superfluid above 
do not respect the law of conservation of angular momentum. Thus, for example, when an 
annular container with superfluid He"^ is cooled, the density of the superfluid is changing; 
and yet it is known that the phase difference (and, consequently, the velocity but not the 
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superfluid angular momentum) remains conserved (Tilley and Tilley, 1986; Donnelly, 1991). 
This is certainly a well-established fact, but it does not settle the basic issue underlying 
the angular momentum generation paradox I have sketched above. In other words, it does 
not say when and to what accuracy the winding number riw along some closed path C 
will be a "good" (conserved) quantum number. For instance, in a bulk superfluid, for a 
fixed closed C, n^, is certainly not conserved in the course of the evolution of the string 
network which follows the quench. This is because strings can move across C, a process 
which obviously alters n^. 

By contrast, in an annulus, cooling (or heating) of the superfluid will lead to a change 
of the superfluid angular momentum (as a result of a change of the Bose condensate 
density ps ~ |*P) but the global conservation law will be nevertheless respected, since 
the superfluid can "push off" the normal fluid, which will in turn dissipate the excess 
momentum on the walls of the container. 

The flrst obvious remark in the wake of the above discussion is to note that, in order 
to be safe from the "bulk" decay of the vortex line network, one must choose the small 
radius of the annulus to be of the order of ^. Otherwise, the resulting "phase around the 
loop" will be modifled by the vortex lines enclosed within the torus. Thus, one will be 
forced to work with capillaries, as the typical values of ^ are close to - 10^ A. This narrow 
annulus may suffer from the possibility of thermally activated transitions (especially near 
Tc), and this process may cause the decay of n^. Perhaps the most likely mechanism for 
such slow decay of the superflow would be due to creation of vortex lines, which can then 
migrate across the annulus, thus changing Uu, by one unit. 

The coincidence of the size of the small diameter of the annulus with the frozen-out 
correlation length ^ contains a crucial clue to the resolution of the angular momentum 
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generation paradox I have described at the beginning of this discussion. This is because 
we can now estimate the angular momentum involved in the superflow as 

Js= (^p.-^7rf -C^ -vr, (4.8) 

where r = C/2tt. The only conceivable source for angular momentum Jg (which, for a 
typical superfluid ps ~ 0.1 g/cm^, | = 10^ A, and C = 1 cm is, in accord with Eqs. (4.4)- 
(4.6), ~ 5 • 10~^^ g cm^/s) would be the angular momentum of Brownian motion, which 
is given by 

Jt= (/esT-M-r^)'/^ . (4.9) 

Here M is the mass of the involved material [which is given by M = • tt^^ ■ C/A in 
Eq. (4.8)]. The above equationcan be readily derived when the kinetic energy of rotation 
(given by J|./(2Mr^)) is equated to the thermal /cbT/2 per degree of freedom. At 
T = T\ = 2°K for the same ps and C we can estimate typical Jy ~ 10~^^ g cm^/s, 
which is much less (several orders of magnitude) than Jg. Thus, at the first sight, it would 
seem unlikely that Jt could be the origin of Jg. 
However, the equation 

Js = Jt (4.10) 

would have to be satisfied not deep in the superfluid phase, where e ~ 1 and ps ~ 
0.1 [g/cm^], as we have assumed above, but, rather, at the freeze-out instant t, for the 
corresponding relative temperature e <^ 1. Now, it is well known that 

Ps = P.(0)e^ (4.11) 

where = 2/3 in the renormalization group theory and in the experiments (Ahlers, 1976), 
and po — 0.36 g/cm.^. In accord with Eqs. (4.8)-(4.10) we get 

pg(7rt) =4{m/hfkBTx, (4.12) 
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which is to be satisfied at e. But the left-hand side of Eq. (4.12) is independent of e as 
^ = ^oc"- Hence, we are led to an unexpected temperature-independent relation: 

2 

TTl 

Psmo = ^^-kBTx/7T (4.13) 

between the density of liquid He^, its correlation length (both at absolute zero), mass of 
the helium atom, fundamental physical constants H and ks, and the critical temperature 
Ta. The left-hand side of Eq. (4.13), for the data cited by Ahlers (1976) is 

Ps{0)Co = 1.44 • 10-^[g cm-2] . (4.14) 

The right-hand side is 

4 (mV^^) ^bT/tt ~ 1.4 • 10-^[g cm-^] . (4.15) 

This coincidence is striking and inspires confidence in the Brownian motion explanation 
of the origin of the angular momentum locked out by the quench. Indeed, Eq. (4.13) has 
been postulated before in the discussion of healing length in superfluid He^ (Ahlers, 1976; 
Ferrel et al, 1969). 

Quench in a Superconducting Loop 

From the experimental point of view, trapping of the winding number by a 
rapid quench in a loop of a superconducting wire may be simpler than the corresponding 
experiment in the superfluid He^. In particular, the detection of the effect — which, in 
a superfluid He*^ involves measurement of small amounts of angular momentum — should 
be now much simpler, as it suffices to measure the magnetic fleld associated with the 
number of ffux quanta trapped during the quench. These flrst impressions may be indeed 
valid as far as the experimental situation is concerned; but from the point of view of 
relating such experiments to cosmology, they are, I believe, too optimistic. As we shall 
see below, quench in a superconducting loop is complicated by the very presence of the 
magnetic (gauge) fleld, which makes the theoretical analysis more difficult and which will 
also complicate experiments. 
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We start by, in effect, repeating the "generic" prediction based on the cosmological 
scenario for the number of fliix quanta generated in converse of a rapid quench: 



= = ^ = (27r)-i ycyl . (4.16) 

Here n$ stands for the number of trapped flux quanta. So far, the effect of the gauge fields 
is ignored. For a loop of radius r = 1 cm and a frozen-out correlation length ^ ~ 10~^ cm 
[see Eq. (3.22)] this yields n$ ~ 3. The corresponding flux would be small but easily 
measurable by the techniques utilizing SQUID's (Tinkham, 1985). 

To arrive at the above prediction, Eq. (4.16), we have paid attention solely to the 
order parameter. But the evolution of the magnetic field trapped by the quench may play 
a role in the final outcome. To consider this possibility, let us first note that the energy 
associated with the trapped fliix $ is nonnegligible: 

= ^^/2L = n| • ^l/2L = nlEo . (4.17) 

Above, Eq stands for the energy associated with a single fliix quantum $o — hc/2e 
[Eq. (2.52)] trapped in a loop with the inductance L: 

Eo = ^l/2L . (4.18) 

In the circular loop of wire the inductance (measured in Henrys) is given by 

L^Att- IQ-^r ln(r/a) [H] , (4.19) 

where r is the radius given in centimeters and 2a is the diameter of the wire. This results 
in typical energies: 

Eo ~ 2.5 ■ 10"^V"^ [erg] , (4.20) 
where we have set r/a = 1000. Energy of thermal excitations is — for comparison — given 

by 

Et ^ ^ksTc ~ 7 X 10-1^ • Te [erg] (4.21) 
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in the vicinity of the superconducting phase transition temperature (which is given in 
°K). 

This comparison — with Et — is relevant because we would like to be able to trace 
the origin of the trapped quanta to either a "cosmological" freeze out of the order parameter 
or to the trapping of the thermal or other fluctuations of the electromagnetic field. This 
second possibility is obviously important, as it is well known that in the course of a 
slow phase transition in an external field the flux trapped by a superconducting loop 
is determined by that magnetic fleld. The above considerations allow us to at least 
differentiate between the experimentally testable consequences of the two alternatives. 
Predictions based on the "cosmological" Eq. (4.16) depend only on the quench rate (which 
determines ^) and on the circumference C of the loop. By contrast, trapping of the thermal 
fluctuations of the field will be affected by the self-inductance L, which can be changed 
without altering C (by, for example, coiling up the loop into a solenoid). Thus, if the 
trapped flux has its origin in thermal fluctuations of the trapped flux, the expected value 
of the fliix will be given by 

By contrast, freeze out of the order parameter leads to the prediction: 

(5$/$o)' = nl = C/i . (4.23) 

The number of trapped quanta should also slowly scale with the quench rate, since n$ in 
the Landau-Ginzburg theory varies as 

n$ ~ Tq^/^ . (4.24) 

The experiment should then be able to tell which of the two alternatives — the one 
driven by the fluctuations of the order parameter or the one in which the gauge field 
dominates — is realized "in practice." This question is of cosmological interest since, as 
we have already indicated, the opinion concerning overwhelming importance of the order 
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parameter in the creation of topological defects is not unanimous (Rudaz and Srivastava, 
1993; Kibble and Vilenkin, 1995). However, it should be noted that there are significant 
differences between the physics which is relevant in the cosmological (bulk) context and the 
one which will play an important role in the superconducting loop experiment considered 
here. In particular, the gauge field in the bulk phase transition acquires mass — and, 
therefore, becomes limited in its range to the penetration depth A = — as the phase 
transition is taking place. Thus, the range over which the order parameter may become 
causally connected will be in any case limited by an exponential to a finite distance, even if 
the phase information could propagate with the gauge field. By contrast, quench in a loop 
leaves the magnetic field inside the loop massless. Therefore, if gauge fields were to play 
a dominant role in setting up relative phases of the broken symmetry Bose condensate, 
the infinite range and the speed with which different sections of the loop can infiuence one 
another could prove significant. 

Detailed consequences of this competition between the order parameter (Bose 
condensate) and the gauge (magnetic) field are difficult to predict in general. We shall 
explore some of the relevant physical phenomena below. We start by noting that, in 
addition to the relaxation time r and the quench time scale tq (which in turn decides the 
freeze-out instant t) , there is at least one more time scale which limits the speed with which 
changes in the magnetic field (or Bose condensate) configuration can take place. The loop 
is an R-L circuit, so the fiux it encloses can change no faster than on a time scale: 

TRL = L/R . (4.25) 

Here R is the resistance. In the vicinity of the superconducting phase transition, this 
time scale is rapidly increasing, as R tends to zero. Moreover, it is dominated by the 
few sections of the loop which are still normal, because of, for example, slightly different 
critical temperature in different sections of the wire. 

This situation (and more generally, any configuration of a two-dimensional ^' in an 
approximately 1-D space) can be represented in a diagram shown in Fig. 9. The quantity 
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of interest in our considerations — known as the fluxoid in the superconductor physics [see 
Eq. (2.55)] — can be used to define the gauge-invariant winding number: 



In effect n^, is the number of times the phase of the order parameter wraps around along 
the loop C. 

In equilibrium, the amplitude of \E' is set by a, Eq. (2.33). The turns are approximately 
equally spaced and can be thought of as regular (although their distribution is gauge- 
dependent). As the phase transition temperature is approached from below, the number 
of turns is constant, but a ~ ^/\e\ becomes smaller, which may eventually allow for changes 
in n^. We shall consider processes which can cause changes of n„ below. Let us begin 
with a situation where the loop C is superconducting except for a small section somewhere 
along its circumference. Then the winding number will not be fixed. 

When a system exists in this state for a time long compared to the time scale on which 
the fiux through the loop can change, the configurations of will explore distribution 
determined, on the one hand, by the energy of these configurations (set both by ^{x) 
and by the magnetic field associated with it) and by the external noise sources which can 
pump energy into the system both through the magnetic field and through When the 
temperature in the normal section of the wire is lowered (so that the two ends of ^{x) are 
joined together), the winding number cannot change any more. 

An experiment which explored this regime (although with a somewhat different 
interpretation in mind) has been actually already carried out. Tate, Cabera, and Felch 
(1984) have heated up a small section of a superconducting niobium loop with a laser. After 
a laser was turned off, the number of trapped fiux quanta inside the loop was measured. 
The random number of quanta had a dispersion corresponding to the temperature 6.78°K, 
which is less than the critical temperature of pure niobium (Tc = 9.17°K). They have 
interpreted this as a result of thermal excitations with the dispersion set by the local 




(4.26) 
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critical temperature in the heated-up section of the loop. The discrepancy in the value of 
Tc was explained as a result of local fluctuation of the critical temperature in the loop. In 
view of our preceding discussion, one might venture an alternative explanation. The lower 
temperature associated with the dispersion of the number of the locked out quanta may 
be a result of the coupling of the system (including the order parameter) with both the 
"hot spot" illuminated by the laser and with the heat bath (which was presumably at a 
still lower temperature than the reported 6.78°K). 
Activation Processes in a Superconducting Loop 

The above discussion has prepared us to consider the problem of the phase transition 
in a loop from the new vantage point with the intuitive understanding based on Fig. 8. The 
first conclusion is already at hand. When the phase transition occurs so nonuniformly that 
the last "break" in the continuity of the Bose condensate disappears only well after the time 
scale needed for equilibration of the rest of the loop (t > rm) we can expect a distribution 
set (or, at least, significantly infiuenced) by the thermal and other fiuctuations driving the 
system. One of these configurations becomes then topologically "trapped" when the whole 
loop becomes superconducting. 

By contrast, when the phase transition occurs simultaneously along the loop, we 
recover the picture we have discussed in the cosmological context. Different sections of the 
loop reach the broken symmetry phase independently. Therefore, at the end of the quench, 
the fiux trapped inside C will be set by the winding number determined by Eqs. (4.1) 
and (4.16). 

It should be noted that the value of the flux is not constrained to the multiple of 
the number of fliix quanta but, rather, by the fliixoid quantization condition, Eqs. (2.55) 
and (4.26). The two are equivalent only in the limit where the penetration depth A is small 
compared with the thickness of the wire, A <^ 2o. When this is not the case (i.e., when 
A > a), the fliix inside C may be quantized in units smaller than $o- This is because in 
that limit both the velocity of the Cooper pairs and the vector potential are important 
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(Blatt, 1961; Bardeen, 1961). This can be investigated in the case where the current is 
approximately independent of the location inside the wire (which is the case when A > 2a). 
The velocity of the Cooper pairs is then given by the gradient of the phase which is not 
compensated by the magnetic flux: 

V = fiq/m* = hq/2m . (4.27) 

Here, m* is the mass of a Cooper pair and q is the residual wave vector: 

Above, the integer n$ is chosen so that qC is in the interval (— tt, tt). 

Equilibrium values of the magnetic flux $ through C will be determined by the minima 
of the total energy associated with a certain flux and with the kinetic energy corresponding 
to the velocity of the charge carriers, Eqs. (4.27) and (4.28): 

E = ^^/2L + Ncm*v^/2 . (4.29) 

Here L is the inductance and Nc is the total number of Cooper pairs in the loop C. It can 
be obtained from the (temperature-dependent) density of the Cooper pairs = I^P- 

Nc = Cira^ • 1*1^ . (4.30) 

Using the above equations, it is straightforward to show that E is minimized when the 
flux is given by 

$ = $0 • n$/(l + Eo/Ek) . (4.31) 

Here, Eq is deflned by Eq. (4.18) and corresponds to a single trapped quantum of flux, 
while Ek is the kinetic energy of the charge carriers for the velocity corresponding to the 
gradient of phase given by 2tt/C: 

^K = iiVc™-(A.|)\ (4.32) 
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It follows that, in equilibrium, the flux in a loop of a given self-inductance L and 
circumference C will be quantized in units of 



$0 = $o/(l + Eo/Ek) . (4.33) 

Furthermore, in the vicinity of Tc the winding number Uw = n$ can be altered by thermal 
activation. The potential barrier which separates different integer values of the winding 
number is approximately given by 

Fb ^ Ek/2 , (4.34) 

for the small values of riy,. This energy can be easily small enough to be of the order of 
UbTc- Thus, in the course of the phase transition and already below Tg, fiy^ can still change 
and may be driven towards a value set by thermal (or other) fluctuations which couple to 
the system. 

The "unwinding" transition we have described above occurs along all of the 
circumference C. By contrast, one can imagine a localized thermal excitation involving 
a correlation-length sized section of the superconducting wire. The free energy barrier 
associated with such a transition is approximately given by the volume in which it takes 
place times the specific free energy of the symmetric state. The exact expression (see 
Chapter 7 of Tinkham, 1985) turns out to be 

AF=^.-.Ae, (4.35) 

where A is the cross section of the wire. Typically AF [Eq. (4.35)] is significantly larger 
than Fbi Eq. (4.34). However, the time scale associated with the thermal excitations within 
one correlation length (given by r) is likely to be much smaller than the time (probably 
of the order of trl) which will set the "whole loop" rate of transitions guarded against by 
the barrier F^, Eq. (4.34), we have discussed before. Thus, it is not clear which of these 
two processes will dominate or whether any of them will have a sufficient rate to reset the 
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number of trapped quanta from the value determined by the frozen-out fluctuations of the 
order parameter, Eq. (4.23), to the estimate of Eq. (4.22), set by thermal fluctuations of 
the flux. 

Last but not least, we note that the study of the dynamics of the order parameter 
in a superconducting loop with one (or more) "gaps" is an intriguing subject in its own 
right. Thus, one could create a number of independent "domains" in the superconductor 
by cutting up the order parameter (i.e., by heating up the loop locally). These 
superconducting domains can be then welded together, and the resulting trapped winding 
number (or the resulting fliix) compared with theoretical expectations, which would have 
to be based on some compromise between the order parameter freeze out and thermal 
fluctuations driven activation process we have discussed above. One can imagine designing 
experiments in which the individual domains are small enough so that Uw is given by 
Eq. (4.16), or alternatively, where the flux is essentially thermal, with the expectation 
value set by Eq. (4.22). It should be perhaps pointed out that such experiments could 
also shed a new light on the fascinating questions about the nature of the phase of the 
wave function of the Bose condensate raised by Anderson (1986) and Leggett (1980) and 
intimately related to the interpretational issues of quantum theory (Wheeler and Zurek, 
1983; Zurek, 1991). 



The main difference between the dynamics of the phase transition in the condensed 
matter and cosmological contexts arises from a different scaling of the relaxation time 
scale with the relative temperature e. Order parameters (fields) which are relevant for 
cosmology obey the equation of motion with the second time derivative on the left-hand 
side. Consequently, relaxation time is 



5. 



COSMOLOGICAL IMPLICATIONS 



1 



1 



(5.1) 



r = 
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(Above, we have returned to^ = c = A;B = l.) One immediate consequence of this scahng 
is that the velocity with which perturbations of the order parameter can propagate, given 
by ^/r, remains finite even at Tc, since 

\/\a\ \/\e\ 



For the purpose of evaluating the initial density of the topological defects in the 
cosmological phase transitions, we should also note that a' = a/e is approximately 

a' ^ PT^ (5.3) 

in field-theoretic models. 

In the early Universe the equation of state is dominated by radiation. Thus, the rate 
at which the phase transition happens will be determined by the relation 

TH = TMpL , (5.4) 

which holds in the radiation-dominated Universe. Above t is the time since the "Big 
Bang," and MpL is the Planck mass. The coefficient T depends on the effective number 
of different spin states s of relativistic particles: 

r = ^-#- (8.5) 

The quench time scale can be then obtained directly from Eq. (5.4) as 

TQ = 1/e = 2r Mpl/T^ . (5.6) 

Moreover, the relaxation time can be expressed as 

T = To/ ^/\e\ , (5.7) 

where tq can be evaluated with the help of Eqs. (5.1) and (5.3); 

To - [VPT.y' . (5.8) 
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We have all of the ingredients required to compute the cosmological freeze-out time £, 
which now obeys the equation 



ro/^t/TQ = t. (5.9) 

This is the relativistic version of Eq. (3.3) we have solved before for superfluids and 
superconductors, but now the solution reads: 

i=r',/%/\ (5.10) 

The difference with the previously obtained i = ^toTq, Eq. (3.4), stems from the different 
scaling of r with e, Eq. (5.1). 

It is now straightforward to evaluate 

i = {2TMpL/Tcf'^ I (f3^'^T^ (5.11) 

and to obtain the corresponding e; 



e = /3-i/3/_^ (5.12) 



ji \ 2/3 

2VMpL ) 

We predict that the freeze out will take place at e and that it will determine the overall 
structure of the initial configuration of cosmological defects. 

The characteristic correlation length would be then of the order of 

The corresponding density of topological defects obtains from the usual argument (i.e., 
one monopole/one | section of a string/one section of membrane per volume of |^). 
This value of the freeze-out relative temperature e is, of course, quite different from what 
obtains when the Ginzburg condition is employed. In that latter case: 

ec ~ . (5.14) 
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This would in turn lead to the correlation length of 



-1 



(5.15) 



which is, again, quite different from the estimate of Eq. (5.13). Indeed, the ratio 



1/3 



(5.16) 



shows that, for the usually assumed small /3, ^ could be smaller than for very high 



temperature phase transitions. However, when Tc <^ MpL — that is, for late phase 



Whether the difference between the resulting initial densities of topological defects 
will be of great significance for a cosmological model is, of course, a separate question. 
For example, in the case of cosmic strings most of the observable consequences will be 
determined by the structure of the string network at late times and will be decided by the 
dynamics of string interactions, etc., rather than by the details of their initial configuration. 
On the other hand, when — as is sometimes proposed — topological defects are used as a 
catalyst in baryogenesis, their density will leave an immediate imprint on the Universe and 
will be critically important. 



We have considered creation of topological defects in course of the second order phase 
transition involving non-conserved order parameter. This selection was dictated by the 
importance of such systems for cosmology as well as by the recent experiment in superfluid 
He^ which afforded a laboratory test of the cosmological scenario. The original cosmological 
motivation notwithstanding, defect-forming dynamics of the second order phase transitions 
is of enormous interest in its own right. 

The key question we have attempted to address concerned the initial density of the 
topological defects. This quantity is a witness to the symmetry breaking dynamics in course 



transitions — ^ is likely to be much larger than $^q. 



6. DISCUSSION 
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of the phase transformation. Two paradigms were put forward to compute initial density 
of topological defects. The older one appeals to the activation process and emphasizes the 
relevance of the Ginzburg temperature — the temperature below which thermally activated 
transitions become prohibitively expensive — on the initial density of the defects. This line 
of reasoning would lead one predict a copious production of defects anywhere between the 
Ginzburg temperature Tq and the phase transition temperature Tg. Moreover, any rapid 
quench trajectory which starts within that region would be expected to lead to more or less 
the same density of defects (set by the correlation length at the Ginzburg temperature). 

In contrast to the thermal activation mechanism sketched out above, one can consider 
a scenario in which the density of defects is determined by the dynamics of the order 
parameter in the immediate vicinity of the critical temperature Tg. The characteristic 
scale (which sets the density of the topological defects) is decided at the instant when the 
critical slowing down makes the order parameter so sluggish that it can no longer keep up 
with the changes of the thermodynamic conditions induced by the quench. This happens 
at the relative temperature e [see Eqs. (3.5) and (5.12) for the typical condensed matter 
and cosmological cases, respectively] . 

Freezeout scenario is guaranteed to set the scale of the initial distribution of 
cosmological defects when e is below the relative Ginzburg temperature eg. This is because; 

|e| > led (6.1) 

implies that the dynamics of the order parameter remains effectively frozen until e is 
reached, so there is simply no time for the activation process. In the laboratory phase 
transition this condition will in principle depend on the quench rate, but is in any case 
likely to be satisfied in the superconductors, and violated in the superfluid He^. In the 
cosmological phase transitions it translates into a relatively simple inequality: 

^2ri;i ■ ^'-'^ 
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This condition follows from Eqs. (5.12), (5.14), and a demand that freezeout should persist 
as the Ginzburg temperature is being traversed during the quench. Thus there is certainly 
a range where the freezeout preempts activation "by default" . 

It may seem more surprising that — as it was anticipated in the original proposal for 
"cosmology in the laboratory" (Zurek, 1984, 1985, 1993), and as it is apparently borne 
out by the experimental results (Hendry et al., 1994, 1995) — the freezeout scenario is valid 
even when inequality (6.1) is grossly violated, so that the activation events are expected to 
occur frequently. Apparently, in the superfiuid (and, presumably, also in the cosmological 
phase transitions) the larger scale set by the critical slowing down is "remembered" , and 
its memory is never erased by the activation occuring on the smaller scales. 

While this may have been a surprise, it was certainly not totally unexpected. For 
one thing, the energetic cost of creating large scale fluctuations is prohibitive. And only 
perturbations which have sizes comparable with the characteristic scale of the string laid 
down during the freezeout can appreciably alter initial density of the string network. 
Moreover, the bulk (~ 70%) of network is expected to belong to a single long string. Its 
geometry may be effected by activation events on the smaller scale corresponding to the 
correlation length at the Ginzburg temperature, but the very fact that the phase transition 
has already taken place implies that such small scale fluctuations do not appreciably re- 
arrange long range order. 

Seen in this light, the apparent lack of importance of the Ginzburg temperature for 
the initial density of cosmological defects does not appear too surprising. The symmetry 
is broken — long range order sets in — at the critical temperature. Topological defects are, 
after all, imperfections of that long range order. Hence, persistence of the large scale 
structure imposed at the time when long range order comes into being for the flrst time 
can be reconciled with the small scale perturbations caused by thermal activation. 

This is not to imply that the above arguments and the apparent accord between the 
experiment and the scenario which appeals to the freezeout (and an even more obvious 
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discord with the consequences of the activation scenario) settles all the issues which are 
worth settling. On the contrary, understanding of the dynamics of the process emergence 
of the long range order in the course of the second order phase transitions is only at its 
inception. Issues which have not been adressed in detail (such as the interplay between 
the freezeout and activation) can be perhaps tackled by numerical simulation (preliminary 
results are already in; Laguna and Zurek, in preparation). They also include mode 
dependence of the critical dynamics (Gill and Rivers, in preparation), the influence of the 
gauge fields, as well as a study of problems specifically designed to address experiments. 
Moreover, experiments themselves should not be restricted to just quenches through the 
critical region. Quenches from just below the critical temperature, as well as transitions 
from the broken symmetry phase to just above Tc followed by a quick return to the broken 
symmetry phase should allow one to probe the timescale on which the memory of the 
broken symmetry is erased. One can also imagine performing quenches into the superfiuid 
starting from within the solid phase of He^ or He^. The order parameter responsible for 
the existence of topological defects (presumably) disappears when superfiuid turns into 
solid. However, second sound in the solid He^ exists, which suggests that at least some 
aspect of the superfiuid behavior may persist in the solid as well. It is of course difficult 
to imagine existence of vortices in the solid. But this (and other similar "out on the limb" 
speculations) can be experimentally addressed by checking that the topological defects can 
be created in course of the solid-superfiuid pressure quench, and by finding out whether 
the topological defects imprinted on the superfiuid can weather the transition into the solid 
phase. (If they did - which seems unlikely - then it should be possible to freeze out and 
then to "defrost" vortex lines.) 

This paper was more a preview of the exciting possibilities rather than a review in the 
traditional sense. Its success should be measured by how much it contributes to its own 
obsolescence. For, it is hoped that the research carried out in the exciting areas outlined 
above will be rapid enough to surpass and eclipse the developments described here. 
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Fig. 1. — Effective potential in Landau-Ginzburg model of second order phase transi- 
tion, V{(fi) = q;|<^P+||<^|^. Above the phase transition temperature Tc (Fig. la) (coefficient 
Q! > 0), there is just a single minimum for V{(p) at {(p) = 0. Below the phase transition the 
minimum is degenerate and corresponds to the expectation value of {(p) = a = y^\a\/P, 
with AV = V(0) - V{a) = a'^/2f3. 

Fig. 2. — Phase diagram of He"^. Liquid Hel is known as "normal," while Hell is 
superfluid. 

Fig. 3. — Structure of the superfluid vortex, (a) The density of the superfluid increases 
outward from a "normal" core on a scale set by the correlation length ^. Velocity of the 
superfluid falls off with the distance so that at every radius the topological constraint 
implied by Eqs. (2.27) and (2.36) is satisfied. This leads to one- dimensional vortex lines, 
shown in Fig. 3b. 

Fig. 4. — Schematic trajectory of the pressure quench which can induce a rapid phase 
transition from normal He^ into the superfluid (a). As the A-line is traversed, equilibrium 
relaxation timescale r, Eq. (2.45) diverges (6). Thus, as the quench proceeds at a flnite 
pace, the order parameter will not be able to adjust any more when the relaxation timescale 
becomes comparable to the time when the critical temperature is reached. This leads to a 
freezeout timescale t; within the time interval [— t, t] the quench is in effect instantaneous. 
Thus, the correlation length ^ corresponding to the instant t will be frozen out and will 
set the density of topological defects. 

Fig. 5. — Formation of topological defects in the course of a rapid phase transition. 
Above the phase transition two-dimensional order parameter fluctuates around the 
minimum of the potential, assuming instantaneously values which "point" in approximately 
the same direction on the {'^x, *j/) plane in the domains of correlation-length size ^. If this 
configuration was frozen out by an instantaneous phase transition, the symmetric vacuum 
would be trapped, resulting in formation of topological defects. 

Fig. 6. — Actual quench trajectories in He^ (after Hendry et al, 1995). Trajectory (A) 
never crosses or even approaches the A-line, and the corresponding quench does not lead to 
vortex line production. Trajectories which cross the A-line [such as (C)] result in copious 
vortex line densities. Trajectories which do not cross the A-line, but come very close to 
it [as is the case for (-B)], also result in vortex line production, perhaps as a result of a 
freezeout of thermal population of vortices which may exist in the vicinity of the phase 
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transition temperature. 

Fig. 7. — Phase diagram of He^ with a schematic indication of a few possible quench 
trajectories. In contrast to He^, normal-superfluid quench in Hehum 3 involves increase of 
pressure. One could also contemplate quenches from solid to superfluid for both isotopes 
pf Helium. 

Fig. 8. — Rapid quench in an annulus. Domain size ^ will now determine the typical 
velocity of the superflow. When ^ > 2r, the problem is reduced to the evolution of a two- 
dimensional (i.e., complex) order parameter in a one-dimensional space and is illustrated 
in Fig. 9. 

Fig. 9. — Evolution of the order parameter in the course of a quench in a thin annulus 
can be illustrated with a version of the Argand diagram shown above. Real and imaginary 
parts of the order parameter are plotted in the radial and in the vertical direction as a 
function of position along the circumference of the annulus. For a suflBciently thin annulus 
order parameter will depend on only one variable — it will vary only with the location along 
the circumference of the annulus, which is shown as a greu "doughnut". An example of 
a possible instantaneous state is shown with the black line. Above the phase transition 
temperature order parameter would fluctuate about the energetically favored = 0, 
changing its value significantly on a spatial scale given by the correlation length ^. However, 
after the symmetry is broken typical value of j^^l will be set by the minimum of the potential 
at (7 = a/IckI//?, and the transitions through = will become unlikely, as they require 
activation energy. Consequently, the initial configuration left by the quench will smooth 
out and be able to partialy unwind (oppositely oriented twists of the spiral will cancel). 
The leftover winding number (given by the nember of times the black line wraps around) 
will be stabilized and will result in a "permanent" superflow (in He^) or supercurrent (in 
superconductors) . 
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